Effect of the measurement on the decay rate of a quantum system 
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We investigated the electron tunneling out of a quantum dot in the presence of a continuous 
monitoring by a detector. It is shown that the Schrodinger equation for the whole system can be 
reduced to new Bloch-type rate equations describing the time-development of the detector and the 
measured system at once. Using these equations we find that the continuous measurement of the 
unstable system does not affect its exponential decay, exp(— Ft), contrary to expectations based on 
the Quantum Zeno effect . However, the width of the energy distribution of the tunneling electron 
is no more F, but increases due to the decoherence, generated by the detector. 
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It was suggested that an unstable quantum system 
slows down its decay rate under frequent or continuous 
observations 0. This phenomenon, known as the quan- 
tum Zeno effect, is believed to be related to the projection 
postulate in the theory of quantum measurements ||^ . In- 
deed, in the standard example of two-level systems, the 
probability of a quantum transition from an initially oc- 
cupied unstable state is Q{At) = a{At)^. If we assume 
that At is the measurement time, which consists in pro- 
jecting the system onto the initial state, then after N 
successive measurements the probability of finding the 
unstable system in its initial state, at time t = NAt, 
is P{t) = [1 - a(Ai)2](*/^*). It follows from this result 
that P{t) 1 for At 0, i.e. suppression of quantum 
transition. 

Originally quantum Zeno effect has been considered as 
a slowing down of the decay rate [|| of quantum systems 
in which a discrete initial state is coupled to a contin- 
uum of final states. This coupling leads to an irreversible 
exponential decay from the discrete state to the contin- 
uum of states. This situation is very often encountered 
in physics, as for instance, the a-decay of a nucleus, the 
spontaneous emission of a photon by an excited atom, 
the photoelectric effect, and so on. But from the theo- 
retical and experimental point of view, the effort has been 
mainly concentrated on quantum transitions between iso- 
lated levels Q characterized by an oscillatory behavior 
between the different states. In this latter case the slow- 
ing down of the transition rate has, indeed, been found. 
However, this was attributed to the decoherence gener- 
ated by the detector without an explicit involvement of 
the projection postulate [§,1). On the other hand, the 
slowing down of the exponential decay rate still remains 
a controversial issue, despite the fact that it is exten- 
sively studied and further investigations are clearly 
desirable. 

In this letter we focus our attention on the quantum 



Zeno effect in exponentially decaying systems, using a 
microscopic description which includes the measurement 
devices. The latter is an essential point missed in many 
studies. This work is motivated by the distinct differ- 
ence between continuum energy levels and discrete lev- 
els, stated above, as well as by the theoretical and exper- 
imental importance of the subject. We also propose an 
experimental set up which is within reach of nowadays 
experimental techniques and within which the quantum 
Zeno effect for exponentially decaying systems can be in- 
vestigated. Our results showed that while the decay rate 
of the quantum unstable system is unaffected by the mea- 
surement the energy distribution of the emitted particles 
can be strongly affected. 
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Fig. 1: The point -contact detector near the quantum 
dot. The energy level Eq of the dot is occupied by an 
electron, which tunnels to continuous states Ea of the 
reservoir, fi^ and fi^j are the Fermi levels in the emitter 
and the collector, respectively. 

Fet us consider an electron tunneling out of a quan- 
tum dot to a reservoir of very dense (continuum) states. 
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Ea ■ The dot is placed near a quantum point-contact con- 
nected with two separate reservoirs (Fig. 1). The reser- 
voirs are filled up to the Fermi levels /i^ and /i^j, respec- 
tively. Therefore the current / = e'^TV/2Ti flows from 
the left (emitter) to the right reservoir (collector), where 
T is the transmission coefficient of the point-contact and 
eV = (hl — tJ^R) is the bias voltage [Q. (We consider 
the case of zero temperature). However, when the dot 
is occupied, Fig. 1, the transmission coefficient of the 
point-contact decreases (T' < T) due to Coulomb repul- 
sion generated by the electron inside the dot. Respec- 
tively, the current through the quantum dot diminishes, 
/' < /. Thus, the point-contact does represent a detec- 
tor, which monitors the occupation of the quantum dot. 
Actually, such a point-contact detector has been success- 
fully used in different experiments [0. Notice that the 
current variation (/ — /') can be a macroscopic quantity 
if the applied voltage V is large enough. The dynam- 
ics of the entire system is determined by the many-body 
time-dependent Schrodinger equation i\'^{t)) — H\^[t)), 
where the total Hamiltonian consists of three components 
H = Hqo + Hpc + H^nt^ describing the quantum dot, 
the point-contact detector, and their mutual interaction, 
respectively. These three parts can be written in the form 
of tunneling Hamiltonians, as 

Hqd = EqcIcq + ^ Eacj^Ca + '^[ilacl^CQ + H.C.], 

Q a 

Hpc = Eicjci + Erclcr + ^[r2;,.c|c,r + H.C.] 
Hint = ^[(5il/rc|cr + H.c]clco (1) 

l,r 

Where the operators cf{ci) correspond to the creation 
(annihilation) of an electron in state i. The and ri;^ 
are the hopping amplitudes between the states Eq, E^ 
and El, Er respectively. These amplitudes are found to 
be directly related to the tunneling rate of the electron 
out of the quantum dot (F) and to the penetration co- 
efficient of the point-contact (T) as F = 2Tr\^la\'^p and 
T = {2Tr)^\Qir\^pLPR, respectively. Here p are the den- 
sity of states in the corresponding reservoirs. The quan- 
tity SQir — ^[j. — ^lir represents the variation of the point- 
contact hopping amplitude, when the dot is occupied. In 
our derivations we assume that f2 and p are weakly energy 
dependent and {hl — P^r) ^ f^^p. The latter condition 
is necessary for the exact solubility of the model. 

Consider the entire system in the initial condition, cor- 
responding to occupied quantum dot and filled reser- 
voirs up to Fermi levels p,L and pp, Fig. 1, denoted 
by |*(0)) = |0). This state is not stable: the Hamil- 
tonian (|l|) requires it to decay to continuum states hav- 
ing the form c^^c\^ ■ ■ • CjCj-^ ■ ■ ■ |0). In general, the total 
wave-function at time t can be written as 



l*(0> 



cLco 



l,r 



balr{t)clclciCQ 



|0>. (2) 



Where b{t) are the probability amplitudes of finding the 
system in the state defined by the corresponding creation 
and annihilation operators. Using these amplitudes one 
finds the reduced density-matrix by tracing out the irrel- 
evant degrees of freedom. This density matrix will give 
us all probability distributions describing the behavior of 
the entire system. For instance, aj^ (t) = |6o(i)P, is the 
probability of finding the system in the initial state at 
time t, cjQQ^t) = X); r \bir{t)\'^ is the probability of find- 
ing one electron in the collector and the quantum dot 
is occupied, aaa{t) = X]q I^q(*)P is the probability for 
the electron to tunnel out of the dot into level Ea and 
no electron arriving the collector, and so on. In general, 
the total probability for the electron to occupy the dot 
is (Too (t) = J2n '^00 (0 ; the probability of tunneling 
into level Ea is (7aa{t) = ^„ credit). Here the subscript 
n denotes the number of electrons reaching the collector 
at time t. The corresponding off-diagonal density-matrix 
elements aoa = cr*Q describe the electron in the linear 
superposition of the states Eq and E^ ■ 

In order to find the amplitudes b(t), we substitute 
Eq. (H) into the time dependent Scrodinger equation and 
use the Laplace transform b{E) = b{t) exp{iEt)dt. 
Then we find an infinite set of algebraic equations for 
the amplitudes b{E), given by 



Ebo - Y ^<yba - Y, ^'ir^lr = « > 
a l^r 

{E + Eq — Ea)ba — f^a^O — ^ ^Irbalr — , 

l,r 

{E + El — Er)bir — ^ ^abalr 
a 

^nlbo-Y^l'r~b^l'rr' =0 



(3a) 
(3b) 

(3c) 



l',r' 



It is very important that the tracing of the reservoir vari- 
ables can be carried out directly in Eqs. (^ without their 
explicit solutions. As a result, Eqs. (||) are converted to 
the Bloch-type equations for the reduced density-matrix 
cr|" -*(/:). Such a technique has been derived in In 
this paper we generalize it by converting Eqs. (^ into rate 
equations without tracing over all the continuum states. 
We, thus, obtain generalized Bloch-type equations which 
determine the energy distribution of the tunneling parti- 
cles. In the following we outline this derivation, relegat- 
ing the technical details to a more extended publication. 
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First, we replace each of the sums in Eqs. by an 
integral, J^k ~^ I p{Ek)dEk which can be treated ana- 
lytically. The Eq. (^a|), for instance, after replacing in it 
the amplitudes ha and hir by the corresponding expres- 
sions obtained from Eqs. (|^), (|3cl), becomes 



E 



E + Eq ~ Ea 
Vl'lpL{El)pR{Er)dEldEr 

E + Ei-Ei-Er 



bo{E)+T = 



(4) 



where J- denotes the terms in which the amplitudes b 
cannot be factorized out the integrals. These terms van- 
ish when the integration limits are extended to infinity 
(the large bias hmit, (/xl — Hr) ^ ^\i^pl,r)- This is 
due to the fact that the singularities of the amplitudes h, 
as functions of the variables Ei^r,a, he on the same side 
of the integration contour. The remaining integrals in 
Eq. (^) can be splited into a sum of singular and princi- 
pal parts. The singular parts yield ir/2 + iD' /2, where 
D' = {2'k)\^1[^\^ plPr{pi — Pr)- While the principal parts 
induce a shift of energy which is merely absorbed by a 
redefinition of the energy levels. Performing the same 
procedure with all other equations we reduce them 
to the following system of equations for the amplitudes 
b{E) |,|| 



D' r' 

E + i—+i-] bo{E)^i , 



D 



(5a) 



E + Eo ^ Ea + I- j baiE) - n^boiE) = , (5b) 
E + Ei-Er + iY+ 4 ) bir{E) - nlbir{E) = (5c) 



where D = {2Tr)\flir\'^ pLPRi^J■l ~ Pr)- In order to trans- 
form Eq. (|^) to equations for the density-matrix we multi- 
ply each of them by the corresponding complex conjugate 
amplitude b*{E'). For instance by multiplying Eq. ( ^ ) 
by b*^{E') and subtracting its complex conjugated equa- 
tion multiplied by ba{E) we obtain 



{E-E' + iD)ba{E)bl{E') 

^nME)bl{E')-bl{E')ba{E') 



(6) 



It is quite easy to see that the inverse Laplace transform 
turns this equation to the following one for the density- 
matrix 



(7) 



Proceeding in the same way with all other equations 
(||) and integrating over the continuum states of the col- 
lector and the emitter, we obtain the following infinite 
set of equations for the density matrix a(i) 



v(") 



-(r + D>oo 



00 



(8a) 
(8b) 



(«) 
00 



T + D + D' 



J") 

^aO 



(8c) 



Eqs. (H) are a generalization of the previously de- 
rived Bloch-type rate-equations for quantum transport 
in mesoscopic systems PJl2|. They have a clear phys- 
ical interpretation. Consider for example Eq. ( ^a[ ) for 
the probability of finding the electron inside the dot and 
n electrons in the collector. It decreases due to one- 
electron hopping to the collector (with rate D'), or due 
to the electron tunneling out of the dot (with rate F). 
These processes are described by the first ("loss") term 
in Eq. (pa|). On the other hand, there exists the opposite 
("gain") process when the state with (n — 1) electrons 
in the collector converts into the state with n electrons 
in the collector. It also takes place due to penetration 
of one electron through the point-contact with the same 
rate D' (the second term in Eq. (^a|)). 

The evolution of the off-diagonal density-matrix ele- 
ments is given by Eq. (^). It can be interpreted in the 
same way as the rate equation for the diagonal terms. 
Notice, however, the difference between the "loss" and 
the "gain" terms. The latter can appear only due to co- 
herent transition of the whole linear superposition [^,0 . 

Since our rate equations distinguish between different 
continuum states (Ea), we can find the energy distribu- 
tion of the tunneling electron by tracing out the detector 
states n in Eqs. (|8|). As a result we obtain the following 
final equations for the electron density-matrix a(t) 

(Too = -FtToo (9a) 

&aa — i^ai^TOa — CTao) (9b) 

F ^ F(/ 

&aO — i{Ea - Ea)(JaQ — if^aCToQ <J aO ■ (9c) 



Here Td — {vD~ vD') is the decoherence rate [g. 

It is instructive to compare Eqs. (^ with the similar 
Bloch-type equations describing quantum transitions be- 
tween two isolated levels fslJl^. In the case of isolated 
levels {El and E2), the equations for the density-matrix 
a are symmetric with respect to Ei and E2. Whereas in 
the case of transition between the isolated [Eq) and the 
continuum states {Ea) the corresponding symmetry, be- 
tween Eq and -Eq, is broken as can be seen, for example, 
in the equation for the off-diagonal term a^a where the 
coupling with aaa is missing. Eq. (p^). 

The probability of finding the electron inside the dot 
is obtained directly from Eq. (pa|) given by 



0-00(0 = exp (-Ft). 



(10) 



It means that the continuous monitoring of the unstable 
system does not slow down its exponential decay. Nev- 
ertheless, it can be shown that the energy distribution of 
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the tunneling electron, P{Ea) = (Taa{t oo), is affected. 
Indeed, by solving Eqs. (j^) in the limit of i — s- cx) we find 
a Lorentzian distribution centered about — Eq'- 



P{Ea.) 



iEo-E^)^ + ir + rar/4- 



(11) 



If there is no coupling with the detector, = 0, the 
Lorentzian width (the line-width) T is exactly the inverse 
life-time of the quasi-stationary state, Eq. (^. However, 
it follows from Eq. (|l^) that the measurement results in 
a broadening of the line-width, which becomes F -I- 
due to the decoherence generated by the detector. At 
first sight this result might look very surprising. Indeed, 
it is commonly accepted that the line-width does corre- 
spond to the life-time. Yet, we demonstrated here that it 
might not be the case when the system interacts with an 
environment (the detector). To understand this result, 
one might think of the following argument. Due to the 
measurement, the energy level Eq suffers an additional 
broadening of the order of F^. However, this broadening 
does not affect the decay rate of the electron F, since the 
exact value of Eq relative to Ea is irrelevant to the decay 
process. In contrast, the probability distribution P{Ea) 
is affected because it does depend on the position of Eq 
relative to Ea as it can be seen in Eq. (|Tl]). 

Although our result has been proved for a specific de- 
tector, we expect it to be valid for the general case, pro- 
vided that the density of states p and the transition am- 
plitude vary slowly with energy. This condition is suf- 
ficient to obtain a pure exponential decay of the state 
Eq. On the contrary, if Pa{Ea) or U.^ depend sharply 
on energy, then the integrals in Eq. (|^) yield additional 
i?-dependent terms that modify both the exponential de- 
pendence of the decay probability, Eq. ( p^ ) and the en- 
ergy distribution, Eq. (|ll| ) . As a result, the measurement 
process could hinder the decay rate |14| . 

We emphasize that our results were obtained from the 
Schrodinger equation describing the dynamical evolution 
of the entire system, without explicit use of the projection 
postulate. This is in contrast with other works, as for in- 
stance where the reduction was repeatedly involved 
during the continuous measurement process. Although 
our final result does not display any slowing down of the 
decay rate, it should not be considered as a contradiction 
with the projection postulate. Indeed, the hindering of 
the decay rate, generated by the projection postulate, re- 
lies on the assumption that the probability of transitions 
between different quantum states is Q(At) = a(At)^. 
It is definitely correct for transitions between isolated 
states, where the transition probability has an oscilla- 
tory behavior. However, in the case of transitions from 
isolated to very dense states, Q{t) represents a sum of 
many oscillations with close frequencies. When averaged 
over the time- interval 5t = '^/£, where £ is the width 
of the function \Vl{E)\'^ p{E) (in our case £ oo), the 
resulting Q{t) would then represent a pure exponential 



decay. In this case Q{At) oc At, so that repeated ap- 
plications of the projection postulate would not change 
the life-time of the decayed state. But, if 5t is finite that 
zero, due to the energy dependence of p or fl, then (Too(i) 
exhibits deviations from the exponential behavior, which 
could result in quantum Zeno effect Q|. 

In conclusion, we have given a microscopic descrip- 
tion of quantum Zeno effect in pure exponentially decay- 
ing quantum systems including the measurement devices. 
Our results show that while the measurements does not 
affect the decay rate, the energy distribution of the tun- 
neling electron is broadened. This description applies to 
a wide range of physical processes as mentioned at the 
beginning of this letter. In particular, it can be verified 
in experiments with mesoscopic quantum dots, by using 
the point-contact detector or an alternative set up. 
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